Gauge-Invariant Quantum Gravity Corrections to Gauge Couplings 
via Vilkovisky-DeWitt Method and Gravity Assisted Gauge Unification 
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Gravity is the weakest force in nature, and the gravitational interactions with all standard model 
(SM) particles can be well described by perturbative expansions of the Einstein-Hilbert action as an 
effective theory, all the way up to energies below the fundamental Planck scale. We use Vilkovisky- 
DeWitt method to derive the first gauge-invariant nonzero gravitational power-law corrections to 
the running gauge couplings, which make both Abel and non-Abel gauge interactions asymptotically 
free. We further demonstrate that the graviton-induced universal power-law runnings always assist 
the three SM gauge forces to reach unification at the Planck scale, irrespective of the detail of 
logarithmic corrections. We also compute the power-law corrections to the SM Higgs sector and 
derive modified triviality bound on the Higgs boson mass. 
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1. INTRODUCTION 

Although gravity, as the weakest force in nature, is 
more perturbabtive than the other three fundamental 
forces all the way up to energies below the Planck scale, 
it was found to be non-renormalizable in the conventional 
sense [1]. But this does not prevent the enormous range 
of successful physical and astrophysical applications of 
the Einstein general relativity of gravitation. In fact, all 
nature's four fundamental forces can be well described 
by the modern formulation of effective field theories [2] , 
with no exception to gravitation [3] . The leading terms 
in the Einstein-Hilbert action, 



S 



EH 



i (i?-2A ), 



(1) 



are just the least suppressed operators in the effective 
theory of general relativity under perturbative low energy 
expansion, where « 2 = 167rG = 16ir/Mp is fixed by 
the Newton constant G (or Planck mass Mp ~ 1.2 x 
10 19 GeV) and Ao denotes the cosmological constant. 

All standard model (SM) particles must join gravita- 
tional interaction with their couplings controlled by the 
universal Newton constant G. It is thus important to 
understand, under the effective theory formulation, how 
gravity corrects the SM observablcs, in connection to 
the other three gauge forces in nature. Robinson and 
Wilczek [4] initiated a very interesting study of gravita- 
tional corrections to running gauge couplings, but it was 
then realized that their calculation by using conventional 
background field method (BFM) [5] is generally gauge- 
dependent and the net result vanishes [6, 7]. 

However, it is important to note that Vilkovisky and 
DeWitt [8] proposed a new approach over the conven- 
tional BFM, especially powerful for analyses involving 
gravitation, which is guaranteed to be gauge-invariant, 
independent of the choices of both gauge-condition and 
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gauge-parameter [9, 10]. The Vilkovisky-DeWitt method 
was recently applied by Toms to study logarithmic cor- 
rections of graviton to the running coupling of QED with 
a nonzero cosmological constant [11] under dimensional 
regularization, and to scalar mass [12]. 

The purpose of the present work is to apply the fully 
gauge-invariant Vilkovisky-DeWitt method [8] for study- 
ing the gravitational corrections to the power-law run- 
ning of Abel and non-Abel gauge couplings. We derive 
the first gauge-invariant nonzero power-law correction, 
which is asymptotically free, in support of what Robin- 
son and Wilczek hoped. We also extend this approach 
for studying power-law corrections to the SM Higgs sec- 
tor and derive modified triviality bound on the Higgs 
boson mass [13], as will be summarized in the last part 
of this paper. The power-law running originates from 
the quadratical divergences associated with graviton loop 
with overall couplings proportional to k 2 . The gravi- 
tational coupling k 2 has negative mass-dimension equal 
— 2 , so the graviton induced loop contributions can gen- 
erate generic dimensionless power-law corrections to a 
given gauge coupling g i , of the form <7 ? -k 2 A 2 , where A 
is the ultraviolet (UV) momentum cutoff. After renor- 
malization one can deduce the generic form of one-loop 
Callan-Symanzik f3 function by general dimensional anal- 
ysis, 



(47T) 



(4tt) : 



(2) 



where \i is the renormalization scale and the coefficient 
a has to be determined by explicit, gauge- invariant com- 
putation of graviton radiative corrections. There is no 
reason a priori to expect a be exactly zero (as stressed 
by Robinson and Wilczek). Our key point here is to 
extract the physical power-law corrections via a fully 
gauge-invariant method a la Vilkovisky-DeWitt [8] . The 
physical meaning of the quadratical divergences in non- 



2 



renormalizable theories was clarified in depth by Velt- 
man [14] and he advocated to use dimensional reduction 
(DRED) method [15] (rather than dimensional regular- 
ization (DREG)) to consistently regularize quadratical 
divergences as d = 2 poles for the Higgs mass correc- 
tions in the SM. Then, Einhorn and Jones made further 
insight [16] that any regularization procedure which pre- 
serves the right number of spin degrees of freedom for 
each field should give the correct results of quadratical 
divergence. This includes DRED but excludes DREG, 
as DREG miscounts the physical spin degrees of free- 
dom for dealing with quadratical divergences [14, 16]. For 
quadratically divergent integrals, except to know that a 
consistent regularization such as DRED exists for them, 
there is no need to explicitly work out these integrals until 
after we finish computing and summing up all their co- 
efficients via gauge-invariant formulation. Then we can 
identify the remaining single divergent integral and re- 
regularize it at d = 4 by placing a common physical mo- 
mentum cutoff; the renormalization will be carried out 
to extract the power-law corrections. (This procedure 
was applied to extract the gauge-invariant quadratical 
divergence in the Higgs boson mass and was shown to be 
regularization-indcpendcnt [16].) We have explicitly used 
DRED method for our analysis (a la Veltman [14]) and 
checked all possible consistencies. 




FIG. 1: Coordinates difference in field space as a manifold. 



variant, and thus T[<p] may depend on the choice of gauge 
condition. 

The way out of this trouble is to replace the coordi- 
nate difference tp 1 — tp 1 by a covariant vector a 1 [tp, 0\ as 
illustrated in Fig. 1, and introduce a connection T l - k in 
field space, we can expand a 1 as, 

a% tp] = (p 4 -if*- \v) k {tp - pY(tp -p) k + ---. (4) 

Thus, the Vilkovisky-DeWitt effective action Tq can be 
constructed as a scalar in field space [8] , 



e ir °^ = ldpLi[p]e*vi[S[p] + C- ll [0\T G M° 3 &^]\ 

(5) 

where Tq j = STq/Sp and a subscript "comma" will 
be always used to denote the functional derivative. The 
coefficient C~ u can also be expanded perturbatively, 



2. GAUGE-INVARIANT VILKOVISKY-DEWITT 
EFFECTIVE ACTION 

The Vilkovisky-DeWitt approach [8] modifies the con- 
ventional BFM in order to build a manifestly gauge in- 
variant effective action. It is noted that a gauge transfor- 
mation corresponds to a field-rcparamctrization in field 
space tp i — > p[ (where we use the condensed notation 
of DcWitt [17], with the subscript i denoting all internal 
and Lorentz indices besides the spacetime coordinates). 
A variation of the gauge-fixing condition for a gauge the- 
ory is equivalent to a change in the external source term. 
This means that although the classical action S[ip] is in- 
variant under field-reparametrization, the effective action 
T[0\ in the conventional BFM is not a scalar. The defi- 
nition of conventional T[tp] , 



dtpn[<p] exp i 



S[<p] + 



Sip 1 



(3) 



includes the external source 5T[0]/6tp l = —Ji, where tp 1 
denotes the background of tp 1 , and dtp [i[tp] is the mea- 
sure of functional integral. For a gauge theory, fi[tp] 
contains gauge-fixing condition and the corresponding 
DeWitt-Faddeev-Popov determinant. From the geomet- 
rical viewpoint with the field configuration space as a 
manifold [8] , it is clear that the difference of two distinct 
points tp 1 — tp % in (3) is not field-rcparamctrization co- 



Cj u [p, tp) = 5)+-W mnj [tp]a m [tp, tp]a n [tp, tp]+- ■■ , (6) 

where R l mn j is the curvature tensor associated with con- 
nection Tj k . Since the i?* ranj -term (or other higher order 
terms in the expansion) already contains two covariant 
vectors a m and <r n (or more), it could contribute to the 
effective action (5) a term which is at least cubic in a 1 , 
and thus will be irrelevant to the one-loop effective ac- 
tion. 

Under perturbation expansion, we can write down the 
one-loop Vilkovisky-DeWitt effective action, which is also 
a scalar under reparamctrization, 



r G [<p) = S[0] -i\nfi[0\ + -TrlnV m V„S, (7) 

where V TO is the covariant derivative associated with 
connection TL„ , 



S 2 S 



5tp m 5p n 



^ s =s^~ T ^W- (8) 



For a gauge theory, consider the infinitesimal gauge 
transformation, 



5tp l 



K' 



(9) 



with K l a [tp] being the generators of gauge transformation 
and e a the infinitesimal gauge-group parameters. Thus, 
for quantization the gauge-fixing condition x a (<p) and 
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the DeWitt-Faddeev-Popov ghost term Q a p[0\ = jjr 
should be introduced. So, the Vilkovisky-DcWitt effec- 
tive action is given by, up to one-loop order, 



S[ip] - Tr In Q a/ 3 [<£] 

+ ^Trln(v m V„S+^( Xct X Q ),^ 



(10) 



which is proven to be invariant under the change of gauge 
condition \ a and gauge-fixing parameter £ [8, 9]. We 
have defined the background field (p l and fluctuating 
field ((? via ip l = (p % + (p l . For calculations in a spe- 
cific gauge theory, the connection r* fc is very compli- 
cated and non-local. But it can be shown that, when the 
Landau-DeWitt gauge condition 



Xa[<P,<f] 



Kai[<Pffl 







(11) 



is chosen, the relevant parts of connection coefficients are 
simply given by the Christoffcl symbol associated with a 
metric G?y in field space [18], 



r' 



-G \Lrtj,k 



G 



G 



jk.e) 



(12) 



In summary, Vilkovisky-DcWitt effective action pro- 
vides a fully gauge-invariant description of the off-shell 
gauge field theories, which is guaranteed to be inde- 
pendent of choices of both gauge-condition and gauge- 
parameter. In the following, we will apply this method 
to analyze the quantum gravity coupled to the Abel and 
non-Abel gauge theories, as well as the SM Higgs sector. 



G 



A^(x)A„(y) 



-gg^6(x-y). 



(16) 



Thus, the relevant part of connection can be derived from 
Eq. (12) under the Landau-DeWitt gauge condition (11), 



x 8{x— y)5(x— z), 



(17) 



r 



A„(a:) 



k(V- 5^-S^ a )S(x-y)S(x-z), 



where we have introduced the symmctrization notation, 
S^Sp = \ (6ffi + Sffi) , and so on. The Landau- 
DeWitt gauge condition should be determined by the 
gauge transformations, 

Sg^ v = -g^ a 9ue a - g va d^ - ^ a 9 a g^ , (18a) 

SAy, = -d„e - A v d^ v - e'd^ , (18b) 

with e a being the infinitesimal parameter of gravita- 
tional gauge transformation and e the infinitesimal pa- 
rameter of U(l) gauge transformation. Then, the gauge- 
fixing functions for photon and graviton fields are 



X 



X u = \ d fya - o d n h 



2° ^ 



(19a) 
(19b) 



3. GRAVITATIONAL CORRECTIONS TO ABEL 
AND NON-ABEL /3 FUNCTIONS 

We start from the classical action of Einstein-Maxwell 
theory, which consists of the Einstcin-Hilbert action (1) 
and 



S 



EM 



(13) 



Since the Vilkovisky-DeWitt method does not require on- 
shell background, we expand the metric around the 
Minkowski background rj^ , 



We further split the gauge field A^ as 
A,, = A,, + a,, , 



(14) 



(15) 



with A^ the background field and the quantum fluc- 
tuating field. As shown in Eq. (12), the connection T l - k 
is determined from the metric Cry defined in the field 
manifold. There is a natural choice [8] of the field-space 
metric G\j with its nonzero components given by 



G 



-g 



2k 2 



{g^g^+g^g va -g^g^)8{x-yl 



where F M „ = d^A v — d y A^ and an overall factor — 2/k 
is factorized out in (19b) for the convenience of normal- 
ization. So the Lagrangian contains the following gauge- 
fixing terms, 



C 



gf 



1 

2C ; 



1 

2| 



X 



(20) 



where £ ( £ ) is gauge-fixing parameter for photon (gravi- 
ton) field, and will be set to zero at the end of calculation, 
as required by imposing the Landau-DeWitt gauge con- 
dition. 

Then, we consider the connection-induced terms in the 
Lagrangian, 



C = --T 9 



s., 



where 
S 



A 



o 



9it* i<p ~ oT 
The ghost part of the Lagrangian is given by 



k2 1"s^9o/3 hpvhap 



(21) 

(22a) 
(22b) 

(23) 



-n^" F a F a/3 + -F^ F uX 
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where 5\ and S\\ are the changes under gauge trans- 
formations with e — rj and e x — rj X , and (77, r/ A ) the 
anti-commuting ghost fields. 

With these we can sum up all required terms for com- 
puting the effective action (10), 



S, 



Q 



S 



EM 



Jd 4 x (C c 



£gf + £ gh ) ■ 



(24) 



The one-loop effective action will be deduced from 
|Trln(S , Q) j i J , according to Eq. (10). This has clear di- 
agrammatic interpretation. For our purpose we are in- 
terested in all the bilinear terms of background photon 
field, which correspond to the one-loop self-energy dia- 
grams listed in Fig. 2. 




(a) \sS (b) 




(c) 



(d) 



(e) 



FIG. 2: Graviton induced radiative corrections to photon self- 
energy, where wavy external (internal) lines stand for back- 
ground (fluctuating) photon fields, double-lines for gravitons, 
dotted lines for photon-ghosts, and circle-dotted line for gravi- 
ton ghosts. 



We note that among all five diagrams in Fig. 2(a)-(e), 
the first two also exist in the conventional BFM or di- 
agrammatic calculation (though the present couplings 
differ from the conventional ones), but the last three 
arise solely from the connection-induced contributions in 
the Vilkovisky-DeWitt formulation. We systematically 
compute these diagrams using the Feynman rules from 
Eq. (24), and extract only the quadratically divergent 
parts of loop integrals, 



(a) = (l + 2C) K 2 (;pV,- PuP*)Za. 



(25a) 



(b) = - 2 - C) «Vv- *VV) I2 , (25b) 

(c) = + ^) «VV- P&) ^ - (25c) 



( d ) = -« 2 (P 2 V,-P /1 1>J2<2, 

(e) = 0, 



where the integral 



d d k 1 



(25d) 
(25e) 



(26) 



is quadratically divergent for d = 4 by power-counting 
and gets regularized for d < 2 via DRED with a sin- 



gular pole at d — 2 [14], though we need not to ex- 
plicitly work it out so far [16]. Since the external gauge 
fields in Fig. 2 carry Lorentz indices, we will encounter 
some terms containing k^k v in the integral of loop- 
momentum k , where according to the DRED we sym- 
metrize k^k v — > rf LV k 2 jd with the metric rf LV defined at 
d = 4 and k 2 /d at d — > 2 for identifying the l/(d — 2) 
poles. The physical picture of the DRED is clear [14-16]: 
the Lorentz indices of loop-momenta in the numerator of 
the integral should be separated into those of the met- 
ric rf 1 " or external momenta at d = 4 to preserve the 
right spin degrees of freedom for each field, and then the 
remaining scalar integral over the loop-momenta is regu- 
larized at d dimension (which can always be reduced to 
Eq. (26) for quadratical divergence). Summing up all the 
self-energy contributions (25a)- (25e), we find that all ^ 
poles explicitly cancel, which is a consistency check for 
Vilkovisky-DeWitt method, and we deduce the net result 
in the Landau-DeWitt gauge, 



(a)+(b)+(c)+(d)+(e) 



15 



k2 (p 2 V^~P^Pu) x 2- (27) 



As guaranteed by the Vilkovisky-DeWitt method [8, 9], 
this is a fully gauge-invariant result. Noting that the 
singular pole of the integral (26) at d = 2 just corre- 
sponds to the quadratical divergence of the same integral 
at d = 4 [14], now we are free to re- regularize the integral 
(26) at d — 4 by placing a common physical momentum 
cutoff A , 



d 4 k 1 
J^fk? 



, A 2 

16tT 2 



(28) 



Thus, we can deduce the QED gauge-coupling renormal- 
ization, 



1 



5 2 (M) ,9 2 (A) 



15k 2 (A 2 -// 2 ) 
128ttV 



(29) 



under the minimal subtraction scheme, and the corre- 
sponding renormalization constant, 



Z g = 1- 



15 K 2 (A 2 -// 2 ) 
256tt 2 



(30) 



where 5(A) = Z g g(p) and fi is the renormalization scale. 

From Eq. (29) or Eq. (30), we finally derive the gauge- 
invariant gravitational power-law correction to the QED 
/3-function, 



A/3(ff,/x) 



15 



1287T 



(« M )9: 



(31) 



15 



which is asymptotically free and gives a = 
Eq. (2). There is no logarithmic graviton correction to 
the gauge coupling /3-function in the absence of cos- 
mological constant [7] , since dimensional counting shows 
that gravitational logarithmic corrections could only 
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contribute to dimcnsion-6 operators such as (D^F^) 2 
rather than the standard dimcnsion-4 gauge kinetic term 
(13). With a nonzero cosmological constant A in (1), the 
graviton- induced logarithmic correction can appear [11] 
since Ao has mass-dimension equal 2 and thus the prod- 
uct Aok 2 provides a proper dimensionless parameter 
for one-loop gravitational logarithmic correction. We 
stress that due to the nonzero new result in the above 
Eq. (31), the graviton-induced leading power-law correc- 
tion will eventually dominate gauge coupling running at 
high scales and always drive gauge unification nearby the 
Planck scale, irrespective of the detail of all logarithmic 
corrections. This will be demonstrated in the next sec- 
tion. 

For comparison, we want to clarify how our above anal- 
ysis differs from that of the conventional BFM (or the 
equivalent diagrammatical) approach. The latter corre- 
sponds to setting all the T l - k related connection terms in 
Sec. 2-3 vanish. In consequence, only the diagrams (a)- 
(b) in Fig. 2 survive. Then, if we apply the usual naive 
momentum-cutoff procedure for quadratical divergence, 
we find that Fig. 2(a) and (b) exactly cancel with each 
other, 

(a) = -(b) = ^^^( P \„- Pll p„)l 2 , 

(a) + (b) = 0. (32) 

This also agrees to the null result of the conventional 
diagrammatic calculation in the second paper of [7]. As 
another check, we can apply DRED method to regularize 
Fig. 2(a)-(b) by setting all connection terms vanish. Then 
we find the two diagrams no longer cancel, but their sum 
depends on graviton gauge-parameter £ and thus non- 
physical, 

(a) = 3(K 2 (p 2 r) lxu -p^ Pv )X 2 , 

(b) = -(3 + ()k 2 ( P 2 ■q tlu -p IJi p v )l2 , 

(a) + (b) = (-3 + 2()k 2 (pX„- P ^)1 2 . (33) 

This shows that even for well-defined gauge-invariant 
DRED regularization for quadratical divergence (a la 
Veltman [14]), it is crucial to further use the Vilkovisky- 
DeWitt effective action (including the connection- 
induced new contributions to Fig. 2(a)-(e)) for ensuring 
the full gauge-invariance. This also explains why the 
gauge- invariant nonzero power-law correction (31) was 
not discovered before. 

As one more consistency check of the present cal- 
culation, we have used another way [11] to compute 
the Vilkovisky-DcWitt effective action in the coordinate 
space and with the aid of Wick theorem. Let us define 
the photon and graviton propagators, 

{ a ^i. x ) a Av)) = D^{x,y), (34a) 



(h^ v (x)h a p(y)) =D^ a (j(x,y), (34b) 

with 

D^(x,y) = j-{0-e- ik ^D^{k), (35a) 

D^ a p(x,y) = JjLtLe-W'-tiD^ik). (35b) 

Thus, we derive the effective action for gauge field, 

zT A = (i<? 2 ) - !<<??) (36) 

where Si and S 2 are the action terms containing one and 
two external fields, respectively. Then we compute the 
effective action for Landau-DeWitt gauge by using the 
Cadabra package [19], and deduce the gauge-part, 

7k 2 r 1 
*r A = —2 2 Jd 4 x -F^F^ , (37) 

corresponding to diagrams in Fig. 2(a)-(c), and the sum 
of Eqs. (25a)-(25c). For ghost-part, we derive 

iT gh = (iS 2gh ) - - (S 2 gh ) 

= -k 2 1 2 Jd 4 x l -F^F» v , (38) 

which corresponds to diagrams in Fig. 2(d)- (e), and the 
sum of Eqs. (25d)-(25c). With these we obtain the full 
one-loop effective action, 

r = r A + r gh = frx ^-f, v f^ , (39) 

from which we reproduce the same gauge coupling renor- 
malization as in (29)-(30) and the same power-law cor- 
rection to the (3 function as in (31). 

Next, we discuss the extension of the above analysis to 
non-Abelian gauge theories coupled to Einstein gravity. 
We first note that a non-Abelian gauge theory adds no 
more graviton-induced one-loop self-energy diagram be- 
yond those given in Fig. 2, except the couplings in these 
diagrams may differ from QED. So, let us inspect the pos- 
sible change in each diagram of Fig. 2 for the non-Abelian 
case. 

First, Figs. 2(a) and 2(d) contain pure gravitational in- 
teractions only, so they remain the same for non-Abelian 
theories. Second, Figs. 2(b) and 2(e) do not change too, 
since both graviton and graviton-ghost carry no gauge- 
charge. So both the gauge fields in and outside the loop 
must share the same "color" and thus no extra summa- 
tion of "color" over the loop gauge-field. Third, Fig. 2(c) 
could receive a change due to possible "color" summation 
over the gauge-loop. The relevant changes would come 
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from two places at one-loop level. One is the gauge-fixing 
(19b), which contributes a term of the following form, 



pa jpb a „b 
b pv b po a a a P 



(40) 



Thus, given the two external background gauge-fields 
(from Fp V Fp a ) for Fig. 2(c), there is no more summa- 
tion over the gauge-indices of the fluctuating gauge-field 
in the loop. The other contribution comes from the con- 
nection term, namely the first term in (21), which would 
contribute to Fig. 2(c) via the form, 



^pu^pa a \ a K 



(41) 



This allows a summation over the loop gauge-indices "6" 
and enhances the contribution by an overall factor of the 
number of non-Abelian gauge fields, which equals TV 2 — 1 
for the SU (N) gauge group. But our explicit calculation 
shows that this connection-induced contribution actually 
vanishes for both Abelian and non-Abelian cases. Hence, 
the conclusion is that our graviton induced power-law 
correction (31) is universal for both Abelian and non- 
Abelian gauge theories. This means that the same co- 
efficient a in (2) holds for all gauge couplings. 



4. GRAVITY ASSISTED GAUGE UNIFICATION 

Gauge coupling unification is a beautiful idea that sug- 
gests the three apparently different gauge couplings of the 
SM (as measured at low energies) would converge to a sin- 
gle coupling of the grand unification (GUT) group [20] at 
high scales. The evolutions of gauge couplings from low 
scale to GUT scale are conventionally governed by the 
renormalization group eqtaions (RGEs) with logarithmic 
running [21]. The precision data show that logarithmic 
evolutions of the three gauge couplings do not exactly 
converge for the SM particle spectrum [22], while the con- 
vergence works fine in the minimal supersymmetric stan- 
dard model (MSSM) with one-loop RG running. But 
more precise numerical analyses including two-loop RG 
running reveal that even in MSSM the strong gauge cou- 
pling a 3 does not exactly meet with the other two at the 
GUT scale as its value is smaller than (a 1; a 2 ) by ^3% 
(a 5(7 deviation); so it is necessary to carefully invoke 
the model-dependent one-loop threshold effects [23] . It 
was also argued that the gravity-induced effective higher 
dimensional operators can generate uncertainties larger 
than the usual two-loop effects of MSSM and thus signif- 
icantly alter the gauge unification [24] . 

With the gauge-invariant gravitational power-law cor- 
rections (31), we can resolve running gauge coupling 
OiQi) = g? (/u)/4tt from the RGE (2), 



-cp 



-cp 



47T 



dx 



(42) 



with c : 



(4tt) 2 



atiQj,) = 



15 



8irMp 



So we further deduce, 



e -cpt, + 



boiai(p ) 
4n 



Jpl x 



1 ' 4tt 



(43a) 



(43b) 



where for the estimate in (43b) we have kept in mind 
that fi <^ Alp (such as the choice of fi = Mz be- 
low), and we also expanded the exponential integral for 
H < c-i ~ 1.3M P ~ 2 x 10 19 GcV (which holds for most 
energy regions in Fig. 3). Eq. (43) explicitly shows that 
the evolution of any gauge coupling a^fi) will be ex- 
ponentially suppressed by exp[— c/x 2 ] , which dominates 
the running behavior for high scales above O(10 -2 Mp) . 
Hence, the universal gravitational power-law corrections 
will always drive all gauge couplings to rapidly converge 
to the UV fixed point at high scales and reach unification 
around the Planck scale, irrespective of the detail of their 
logarithmic corrections and initial values. This feature is 
numerically demonstrated in Fig. 3 by using the evolution 
equation (43a). 

In Fig. 3(a)-(b), we have analyzed the gauge cou- 
pling runnings for both the SM and the MSSM, where 
the conventional coefficients &o; in the one-loop RGEs 
are, (6 i, b 02 , b 03 ) = (-%, f , 7) for the SM and 
{boi, &02, &03) = (-f , -1, 3) for the MSSM. (Here 
we have adopted the conventional MSSM particle spec- 
trum without including the superpartner of graviton - 
the spin- 1 gravitino, whose potential contribution to the 
gauge coupling running is worth of a future study.) We 
have also input the initial values, a^ 1 = 59.00 ± 0.02, 
a^ 1 = 29.57±0.02, and aj 1 = 8.50±0.14 , at the Z-polc 
li = M z [22]. 

For the graviton-induced one-loop /3-function (2) with 
a n given in (31), we note that at the Planck scale 

,2\ 



/.t = Mp , the loop-expansion parameter, ^yi (k 2 ^ 2 ) = 
If ]^ - 0.6 < 1, and the condition J^(kV) < 1 

holds for /i < 2 x 10 19 GcV, which is better than the naive 
expectation. Fig. 3 shows that above 10 19 GeV the gauge 
couplings rapidly converge due to the exponential sup- 
pression in (43a); and quite before approaching the UV 
fixed point the three curves become indistinguishable for 
q, < 0.01 (SM) and a, < 0.02 (MSSM), corresponding 
to a scale of about 2 x 10 19 GeV. Given the experimental 
error of a 3 {Mz) = (8.50 ± 0.14) -1 and potential higher 
loop effects, this is enough for a possible consistent unifi- 
cation with a finite value of coupling and at a single scale 
around the Planck energy. For the perturbative one-loop 
running of gauge couplings, the renormalization of k 2 
belongs to higher order effect here and is thus not in- 
cluded. It is useful to note that there are clear evidences 
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FIG. 3: Gravity assisted gauge unifications in the SM 
[plot-(a)] and in the MSSM [plot-(b)]. In each case the 
graviton-induced universal power-law runnings become dom- 
inant above 10 18 GeV and always drive three gauge couplings 
to rapidly converge towards the UV fixed point, nearby the 
Planck scale. 



tal error of a 3 (M z ) = (8.50 ± 0.14)- 1 , this is enough 
for a possible consistent unification with a finite value of 
coupling and at a single scale around the Planck energy. 
Also, taking into account of higher loop corrections and 
possible threshold effects may well fill in any remaining 
tiny gap among aj(/z)'s at this scale. For the MSSM, 
as shown in Fig. 3(b), one needs not to worry about 
the model-dependent threshold effects or the two-loop- 
induced non-convergence around the scale of 10 16 GeV, it 
is quite possible that the GUT does not happen around 
the scale of 10 16 GeV, as in the SM case. Instead, the 
real GUT would be naturally realized around the Planck 
scale, and thus is expected to simultaneously unify with 
the gravity force as well. This also removes the old puz- 
zle on why the conventional GUT scale is about three 
orders of magnitude lower than the fundamental Planck 
scale. Furthermore, the Planck scale unification helps to 
sufficiently postpone nucleon decays, which explains why 
all the experimental data so far support the proton sta- 
bility. In addition, this is also a good news for various 
approaches of dynamical electroweak symmetry break- 
ing [27], such as the technicolor type of theories, which 
often invoke many gauge groups with new strong forces 
at the intermediate scales and makes one worry about 
whether a gauge unification could ever be realized at cer- 
tain high scales in these theories. Fortunately, the uni- 
versal gravitational power-law running for all gauge cou- 
plings found above should drive a final unification at the 
Planck scale. It suggests that the Planck scale unification 
may be a generic feature of all low energy gauge groups 
and is fully consistent with the experimental evidences of 
proton stability. It is also an appearing feature of having 
the possibility of a simultaneous unification of all four 
fundamental interactions around the Planck scale. 



supporting Einstein gravity to be asymptotically safe via 
the existence of nontrivial UV fixed point in its RG flow, 
so it may be UV complete and perturbatively sensible 
even beyond the Planck scale [25, 26]. To be conserva- 
tive, we would consider the one-loop prediction of Fig. 3 
as an instructive extrapolation of the Einstein gravity, 
since this provides an encouraging insight on the impor- 
tant role of gravity for realizing the gauge unification, 
and should strongly inspire and motivate more elaborate 
investigations along this direction. 

With these said, we make a few more comments. 
Fig. 3(a) shows that for the SM gauge coupling evolu- 
tions, despite their familiar non-convergence in the re- 
gion of I0 13-17 GeV, the three couplings do unify around 
the Planck scale. Mathematically, as shown in Eq. (43), 
to make oti(fi) reach the UV fixed point would require 
fi — > oo ; but practically, due to the one-loop exponen- 
tial suppression in (43), the differences among ai(/^)'s 
will rapidly decrease to be smaller than the experimen- 



5. GRAVITATIONAL CORRECTION TO HIGGS 
BOSON COUPLING AND MASS 



Without losing generality, we first consider a real scalar 
field which minimally couples to Einstein gravity, 



(44) 



where V(<j>) is the scalar potential, V = ^m 2 (j) 2 + 4<; 



with 



^ ( m < ) corresponding to the unbro- 



ken (broken) phase. To compute the effective potential 
of scalar field, we expand the graviton and scalar fields 
around their backgrounds, 



9 '[us 



(45) 



We impose the Landau-DeWitt gauge condition for the 
gauge-fixing, 



1 K 



(I, 



(46) 
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Then, we derive the connection-induced terms in the 
graviton-scalar sector of the Lagrangian, 



with 



-KS.thT 



4> fy,„ , 



:S., 



r 9^ 



d P m ~ v(4>) 



-d 2 



6 



(47) 

(48a) 
(48b) 



In (47) the graviton-field connection r s ^" ff was given 
in Eq. (17), and the scalar-field related connections are, 




(a) 



(e) 



(b) 



(c) 



(d) 





(f) 



00 



FIG. 5: Graviton-induced corrections to quartic scalar vertex: 
dashed-line is fluctuating scalar field, double-wavy-line de- 
notes graviton, and circle-dotted-line depicts graviton-ghost. 



0(y)0(2 



-9u.J{x-y)8{x- 



r 



0(x) 

<t>(y)<t>(z 



0(y)g^O) 



0, 



(49a) 
(49b) 
(49c) 



which are derived from the metrics G 



0O)0(y) 



V-g$( x -y) and G g^{*)g a p{y) [as in Eq. (16)]. The 
Feynman diagrams for the graviton-induced scalar self- 
energy corrections and quartic vertex corrections arc 
shown in Fig. 4 and Fig. 5, respectively. 



(a) 




o°o 
o o 



(c) 



(d) 



FIG. 4: Graviton-induced self-energy for scalar field: dashed- 
line is fluctuating scalar field, double-wavy-line denotes gravi- 
ton, and circle-dotted-line depicts graviton-ghost. 



For the scalar self-energy corrections, we compute the 
quadratic divergent part for each diagram in Fig. 4, 



Fig.4(a) 
Fig.4(b) = -3K 2 m 2 I 2 



1 \ , 2(1 M 2 1 2 2 

-A + k \ w Km 

2 V« 4C / 4 



Fig-4(c) = (-1 + ^) 

Fi g .4(d) = K y i 2 , 



X 2 , (50a) 
(50b) 
(50c) 
(50d) 



where the A-term in (50a) comes from the scalar quar- 
tic self-interaction alone, which we have included for the 
comparison with graviton induced corrections and for the 
convenience of analysis. We also note that unlike the case 
of photon self-energy, the external lines of the scalar self- 
energy carry no Lorentz index, so the corresponding loop 
integrals can be directly reduced to a scalar integral with- 
out k^h v like combination of loop- momenta, and thus no 
symmetrization of loop-momenta is needed. 

For the graviton-induced vertex corrections in Fig. 5, 
only the first two diagrams have quadratical divergence, 



Fig.5(a) = -h,n 2 l 2 , 
Fig.5(b) = -3Ak 2 Z 2 , 



(51a) 
(51b) 



and all other diagrams contain logarithmic divergence at 
most. 

Summing up relevant contributions we deduce the two- 
point proper self-energy and four-point proper vertex for 
scalar field, 



r 2 (p) 



\\ + \{ P 2 -26mV 



13 



X 



2 • 



r 4 (p) = — iX - —^-Xk 2 X-x , 



(52a) 
(52b) 



where we do not include the graviton-induced logarithmic 
divergent terms (as in [12]) because they are negligible as 
compared to the dominant power-law corrections. Note 
that the sum of (50a)-(50d) explicitly proves the exact 
cancellation of the \ gauge-parameter poles, which is a 
consistency check of our Landau-DeWitt gauge calcula- 
tion. 

To fully understand the results (52a)- (52b), we have 
compared them with those in the conventional BFM 
(or the equivalent diagrammatical) approach where all 
connection-induced new terms are set to zero. The find- 
ings are summarized in Table- 1. 
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TABLE I: Graviton-induced power-law corrections to scalar self-energy and quartic vertex: Summary of the comparison between 
the conventional approach (with vanishing connection and with general gauge-parameter and Vilkovisky-DeWitt approach 
(VDA) (with nonzero connection, and in Landau-DeWitt gauge with f — ¥ in the end). In each entry of the contribution, a 
common factor k 2 I2 is factorized out. Our summed results of T2 and T4 on the 3rd column agree with Eqs. (52a)-(52b). 



From Eqs. (52a)-(52b), we derive the renormalization 
for scalar coupling, A(/^) = Z 2 Z^ 1 X(A) , with the fol- 
lowing renormalization constants, 

1 



Z A = 1 + J^ 2 ( A 2 -, 2 ), 



(53a) 
(53b) 



we compute the graviton-induced scalar /3-function, 



(54) 



which is not asymptotically free, contrary to the gauge 
coupling /3-function (31) we derived earlier. The pure 
scalar loop correction is logarithmically divergent and its 
renormalization gives the usual non-asymptotically free 



under the minimal subtraction scheme. Then, with (53) scalar /3-function (3q = + 



3X Z 



r 



From the two-point proper self-energy (52a), we fur- 
ther perform the renormalization for scalar mass in the 
on-shcll scheme, which fixes the mass counter term, 



5m 2 



1 



32tt 2 



-X 



25 



2 2 
m k 



A- 



(55) 



where we have defined the renormalized mass, m 2 = 
TOq — Sm 2 with m denoting the bare mass parameter 
in the original Lagrangian, and also included the contri- 
bution from the pure scalar loop. Comparing the two 
terms on the right-hand-side of (55), we note that the 
graviton-induced quadratical divergence is actually much 
softer since the product k 2 A 2 = 16tt(A/M p ) 2 = O(10 2 ) 
for an ultraviolet cutoff A ~ Mp . 

It is straightforward to extend the above analysis to the 
SM Higgs boson, since graviton coupling to scalar fields 
is universal. Let us write down the SM Higgs doublet, 



1 / tt 1 + iir 2 
y/2 \ a + in 3 



(56) 



where = a 2 +^~"^ 7r 2 and a = a + a , with a being 

a=l 

the SM Higgs boson and ir 1 2 3 the would-be Goldstonc 



bosons. The Higgs background field a will equal the 
vacuum expectation value ( VEV) , v , at the minimum of 
Higgs potential. 



Consider the SM Higgs potential, V 



J ($t$) 



A($ f $) 2 , with m 2 < and the VEV, v = v /-m 2 /A. 
Then we recompute the gravitational power-law correc- 
tions to Higgs boson self-energy and quartic vertex in 
Fig. 4-5, and deduce the following, 



l\ + l(p 2 -8m 2 H )n 2 



To 



r ff4 b) = -«A-24Ak 2 I 2 



(57a) 
(57b) 



With these we derive the graviton-induced contributions 
to the Higgs boson /3-function and the Higgs mass counter 
term, 



A/3(A,/i) = +-pi(«V), 



3A_ 

87T 2 ' 
1 



'H 



32tt 2 



T2A + 7m^K 2 ) A 



(58a) 
(58b) 



where we have used m H to denote the physical Higgs 
boson mass. We see that (58a) happens to be the same 
as in (54) while (58b) has different coefficients from (55). 
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As compared to the single real scalar case, the changes 
in the computation of (58a)- (58b) arise from two sources: 
(i) the Feynman rule for the quartic Higgs vertex in the 
SM has an extra factor 6 relative to that from the real 
scalar potential below Eq. (44) ; (ii) there are additional 
Goldstone loops in the SM which contribute to Fig. 4(a) 
and Fig. 5(a). 




10° io 2 10 4 10 6 10 s io 10 10 12 10 



A (GeV) 

FIG. 6: Gravitational power-law corrections to the triviality 
bound on the SM Higgs boson mass m H , where the region 
above each curve is excluded. The red curve includes the 
graviton-induced corrections and the blue curve depicts the 
bound with SM interactions alone. 



where e 2 K m n ~ 1 holds to high accuracy due to the 

2 

tiny factor ^-^m 2 ^ = f tjt — . It is clear that in the 
a — >• limit, the condition (62) reduces to the familiar 
triviality bound in the pure SM [28], 

A 2 8tt 2 v 2 



In 



< 



(63) 



' H 



Keeping this in mind, it is instructive to rewrite our 
graviton-corrected triviality bound (62) as follows, 

- A 2 8ttV 

(64a) 



111 ■ 



< 



X 



' H 



A r 



3(1 + A) 



dx L A 
— / In — 

X m 



> 0. 



(64b) 



H 



where an overall factor e 2 K m n ~ 1 on the RHS of 
(64a) is safely ignored. We find that, because the inte- 
gral X > generally holds, the gravitational power-law 
corrections always reduce the RHS of (64a), and thus 
further tighten the triviality bound relative to (63) of the 
pure SM. The graviton-induced corrections play a dom- 
inant role to enhance the triviality bound for the cut- 
off scale A ~ Mp , as clearly shown in Fig. 6. A sys- 
tematical expansion of the present section (including the 
power-law corrections to Yukawa couplings) will be given 
in Ref. [13]. 



Note that the conventional Higgs /3-function in the 
SM receives only logarithmic corrections, /3o(A, fi) = 
+ 2^2- A 2 , and is not asymptotically free. So, with (58a) 
we can write the summed Higgs /3 function, 



/3{\,fi) = b Q \- + a \n n , 



(59) 



with (b , a ) 



.27T 2 



the running Higgs coupling A(/z) 



g^r) > . Solving (59) we deduce 



■ill*- 2 ,, 2 

e 2 K ^0 



A(m 



2a„2^2 dx 



e 2 



— ■ ( 60 ) 



We see that when the right-hand-side (RHS) of (60) van- 
ishes, the renormalized coupling A(/ii) blows up at the 
Landau pole fi = A^ , 



A l (fi ) = ~ h v e ~ 



-■ (61) 

x 



This means that the SM as an effective theory must have 
an UV cutoff A < A^ . For a given Higgs boson mass 
m 2 H = 2\{m H )v 2 , let us set pi = m H . Thus, from (61) 
we can derive the triviality bound for Higgs boson mass, 

-1 



2v 5ft K 2 2 
< g 2 K m H 



e 2 



K 2 X 2 dX 



(62) 



6. CONCLUSIONS 

The fundamental gravitational force universally cou- 
ples to all the SM particles and can be described by 
the well-defined perturbation expansion in the modern 
effective theory formulation [3] . The Vilkovisky-DcWitt 
method [8] profoundly modifies the conventional BFM, 
and provides the manifestly gauge-invariant effective ac- 
tion for reliably computing quantum gravity effects. In 
this work, we used the Vilkovisky-DeWitt method to de- 
rive the first gauge-invariant nonzero gravitational power- 
law corrections to the running of gauge couplings. We 
found the gravitational power-law corrections to be uni- 
versal, making both Abel and non-Abel gauge couplings 
asymptotically free [cf. Eq. (31) and analyses at the end of 
Sec. 3]. We have demonstrated that the graviton-induced 
power-law runnings always drive the three SM gauge 
forces toward to the UV fixed point, reaching final unifi- 
cation at the Planck scale and irrespective of the detail of 
logarithmic corrections (cf. Fig. 3). This raises the con- 
ventional GUT scale by three orders of magnitude, and 
opens up a natural possibility of simultaneous unification 
of all four fundamental gauge forces at the Planck scale. 
We further analyzed the power-law corrections to the j3- 
function and mass of the SM Higgs boson [cf. Eqs. (58a)- 
(58b)]. We found that the graviton-induced scalar f3- 
function is not asymptotically free, and therefore further 
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tightens the triviality bound on the Higgs boson mass, as 
shown in Eq. (64) and Fig. 6. Further extensions of the 
present analysis for computing the power-law corrections 
will be given elsewhere [13]. 
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